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Abstract 

We present the results of a new perturbation calculation in polymer statistics 

which starts from a ground state that already correctly predicts the long chain 

length behaviour of the mean square end-to-end distance {R%) , namely the 

solution to the 2 dimensional (2D) Edwards model. The {R'ff) thus calculated 

is shown to be convergent in N, the number of steps in the chain, in contrast to 

previous methods which start from the free random walk solution. This allows 

us to calculate a new value for the leading correction-to-scaling exponent A. 

Writing {R%) = AN^^il + BN'^ + CN'^ + ...), where v = 3/4 in 2D, our 

result shows that A = 1/2. This value is also supported by an analysis of 2D 

self-avoiding walks on the continuum. 
36.20.Ey 64.60.Fr 
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I. INTRODUCTION 



The central quantity of interest in polymer statistics is the mean square end-to-end 
distance {R%) for a chain of links. The excluded volume effect causes a 'swelling' of the 
chain beyond the {R%) oc relationship of a free random walk. Thus 

{R%) = AN^''{1 + BN~'^ + CN"^ + ...) , (1) 

in the large N limit, where u is the leading scaling exponent, A,B,C are excluded volume 
dependent coefficients and A is the leading correction-to-scaling exponent [|T|. It is now 
firmly established that in two dimensions (2D) i/ = 3/4 is exact. Despite this, there is 
very little agreement about the value of A. Nienhuis predicts A = 3/2, while Rapaport 
has argued that there is no need for a correction term other than the analytic correction, 
i.e., A = 1. However, many numerical studies have disagreed with these results, with 
estimates for A of 1.2 0, 0.84 p[ and 0.65 0,1,0. These numerical estimates are based 
on results obtained from self-avoiding walks on 2D lattices. Theoretical results are also in 
disagreement. Besides Nienhuis's prediction, which relies on a mapping to an exactly solvable 
solid-on-solid model on the honeycomb lattice. Baker et al [|I0| predict A = 1.18 using RG 
arguments, while Saleur ||ri| predicts A = 11/16 by conformal invariance. Interestingly, 
Saleur also gives evidence for A but he then rejects this result. Perturbation expansion 
techniques [0,|l3, which start from the free random-walk solution, have also been used to 



predict {R%), but these methods have resulted in series which are divergent in and v, 
the excluded volume parameter, and hence a value for A cannot be predicted. The obvious 
confusion in both the numerical and theoretical estimates for A motivates the search for a 
better perturbation procedure as discussed here. 

In this paper we report results based on a new perturbation method, which unlike pre- 
vious studies, starts from a ground state that already correctly predicts the exact large N 
behaviour in 2D, namely the Edwards self-consistent solution |]T^. Our approach mimics 



the cluster integral method in statistical mechanics 15 1. Although it has been shown |16 
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that the Edwards solution cannot be the correct form for the self-avoiding random walk 



end-to-end distribution function [jT^, it does predict an exact exponent z/ in 2D, although 



not in 3D. As far as we know, the reason for the coincidence of v in 2D predicted by the 



self-consistent theory with the recently proved exact value |Tl| is not well understood. Thus 
it appears that there is essential physics in the 2D Edwards solution that is not present in 
3D, as easily seen via a simple dimensional argument (see section IV). Mathematically the 
Edwards solution also has convenient features that we can exploit in a novel perturbation 
expansion. 

In section II, we spell out the method for the approach. Here the essential physics consists 
of (a) perturbation about a new ground state, i.e., about the Edwards solution which is 
known to predict an exact exponent z/ = 3/4 in 2D, and (b) a new perturbation potential 
V{r), the latter being much better behaved than the bare excluded volume potential (see 
Eq. d^)). In section III, we present details of the calculation and results and in section IV 
we conclude with some discussions. 

II. METHOD 

We start from the path integral representation of the exact end-to-end distance distri- 
bution function for the Edwards model |jl8[. This is given by 



/■r(L)=R / 1 /.L f)r( ^\ rL /-i \ 

G{R,L)= D[r]exp --/ ds{-^f - uj ds ds' 5^[r{s) - v{s')]] , (2) 

Jr(0)=O V £ Jo OS Jo Js J 



where L is the total chain length L = Nl, I being the step length of one link, and u = v/l"^, 
where v is the excluded 'volume'. It is known that 

v = e-"('-»^)/'=«^]c/V,, , (3) 

where u{vij) is the pair potential between the zthe and jth segment, is the Boltzmann 
constant and T the temperature. Traditionally two problems arise in dealing with this 
intractable path integral. Firstly, divergences appear in the calculation which must be 



3 



handled carefully, and secondly, the resulting series expansion is a power series of increasing 
L and uj. This leads to a result whose radius of convergence for uj diminishes to zero as 
L increases. This property is the hallmark of modern critical phenomena theory, whose 
resolution was offered by the renormalization group approach ||T^J2^. Historically Edwards 
avoided the divergence problems of such an approach by replacing the point contact potential 
by a self-consistent field W{r) which in 2D is equal to vp{r)/l'^, where p{r) is the one-particle 
potential proportional to r"^/'^ [Q. Therefore 

W{r) = Cv^/^r-^/^ , (4) 

where C = (-\/3/47r)^/^Z~^/'^. Thus the Edwards Green's function Ge(R,L) becomes 

r(L)=R ( I , (dvis)V 



''-e^- ^' = X,„,.o ^i--' [-1 L ''\-ir)- L '^<^'* ) ■ 

Our approach relies on obtaining a better first order perturbation expansion by starting 
from the Edwards ground state and then perturbing this with the difference between the 
self-consistent field and the true point contact potential defined by V{s) below. Thus 

G(R,L) = Ge{R,L) X exp W{s)ds - uj J^^ ds j\s'6'^[r{s) - r{s')]^ 

= GeCR, L) X exp (^j^ V{s) ds^ . 

To a first order perturbation expansion in V{s), this now becomes 



(6) 



G(R,L) = Gij(R,L) X 



1+ V{s)ds + 0{V{s) 
Jo 



= GE(R,L) + Gy(R,L) + ... . (7) 

The Fourier transform of Eq. (j^) is 

G(k, L) = GE(k, L) + Gv{Yi, L) + ... , 

= Gs(k, L) + G^{k, L) + G2(k, L) + ... , (8) 

where G£;(k, L) is the Fourier transform of the Edwards Green function, G'i(k, L) the trans- 
form of the self-consistent field acting on G'£;(k, L), and 62 (k, L) is the transform of the point 



contact potential acting on GeO^, L). If we now define the following Laplace transformation 
function 



dL e-^^Gi5(k, L) = C{Ge{K L)} 



(9) 



then it can be shown [jT3[ that due to the properties of the delta function in Eq. (0) and the 
factorizability of the path integral of Eq. (||) 



G2{k,E) = GE{k,E) 



2 / d"^ 

47^2 



Ge{p,E) 



(10) 



At this stage we note a simplification over the method in ref. []13[ which bypasses the need 
to evaluate the integral of Ge{p,E) as it is non-trivial in our case. This can be seen by 
noting that 

J±^GE{p,E)d^ = |e^P-^G^(R,£;) d'Rd^, 

= JJ^e'^-^d^GE{K,E) d^R, 
= J 6\K) Ge{R,E) d^R , 

= Ge{0,E). (11) 

Thus G2(k, E) is given by 



(12) 



G2ik,E) = GE{k,EYGE{0,E) , 

= GEik,EYC{GEiR,L)n=o} ■ 
Using the convolution property of the Laplace transform, then 

62(1^, L) = [G^(k, EYC{Ge{R, L)n=o}] , 

[GE{k, L - uf] {GeCR, ^x)r=o) du , (13) 

which is a convenient result. Once 6*2 (k, L) is obtained from ( |13|) we can then calculate the 
mean square end-to-end chain distance by 



= J d^RR^G{K,L) I J d^RG{R,L) 

" d 



-4 



dk' 



:Gik,L) 



'G{0,L) . 



(14) 



J k=0 



The formula (|T^) overcomes the need to calculate the integral in Eq. ([TOD which is non- 
trivial given the form of Ge{p, E). Moreover, it is worth noting that the above calculations 
are only tractable in our case because of the essentially (shifted) Gaussian nature of the 2D 
Edwards solution and the trick given in Eq. (p^). It appears that a similar method is not 
possible with more sophisticated ground states |T^. 



III. CALCULATIONS AND RESULTS 



We shall first show that the above method works by following the traditional approach 
and starting from the free walk solution in 2D. Given that 



6*0(1^, L) 



1 / 1 



thus 



£{Go(R,L)r=o} = /: 



XnlL 



and 



Go{k,E) 



Using Eq. (|T2D, Eq. ([TD and Eq. in Eq. ([T|) we obtain 





I 




'2lu; n 




£-1 


' 1 " 
_E 


-£-1 







Introducing a cutoff e at the small L limit, this reduces to 



l-^[Llog(L)-L-Llog(e)] 



IL 1 



27l 



L + 0{iJ 



(15) 



(16) 



(17) 



(18) 



(19) 



which agrees with the previously known result 1T^|T3 ■ 

It can be shown that in 2D the Edwards solution for the probability distribution for 
the end-to-end distance for a chain in the large L limit, i.e., the solution of Eq. (^, is given 
by 
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where 



Ge(R,L) = N{L)exp{-j{R-AL^/y) , (20) 



and A^(Iv) is the normahsation. Unhke Edwards ||Tj], the normahsation is needed in our 
calculation. This is determined by the identity 

J GE{'R,L)d^R=l . (22) 

Given the asymptotic form of GEiRyL) in ( pO] ) we can calculate N{L) from (p^), see 
Appendix A. The result is 

For later reference we note that the exponentially decreasing term in the denominator of 
the above expression for N{L) appears as a natural 'cutoff' and prevents many subsequent 
integrals from diverging. It plays a similar role as the cutoff e, but does not need to be 
artificially introduced as in the perturbation expansion based on the free-walk solution 
given above. 

From Eq. (|lj) it can be seen that only the small k behaviour is required. Thus the 
Fourier transform of Ge(R;L), see Appendix B, becomes 



.2 r 



GeO^,L) = 1-- 
The Laplace transform of the above gives 



2B . 



(24) 



GEik,E) = C{l-^[A^Ly^ + ^L]}, 

By examining Eq. (|^ and the form of the one-particle potential W{r), it is seen that 
the term G'i(R, L) is calculated simply by differentiation of Ge(R,L) w.r.t. the excluded 
volume parameter. Thus 



Gi{k,L) = -v 



2/3 _ 



,2/3. 



d 

9(t;2/3) 
d 



GE{k,L) 



(26) 



9(t;2/3) 

The higher order terms present all decay exponentially and thus are insignificant in the large 
L limit. Since only A has a v dependence, then from Eq. ([2l| ) 

3 



G,{k,L) = --A'L'/^ 



(27) 



We note this is just — |G'£;(k, L) to leading order in L, and it acts to decrease the value 
of {R%). The term G2(k,L) is calculated using Eq. Eq. (|0D and Eq. (H) to order fc^, 



thus 



G2(k,B) 



^2 4 V 2 5 



/:{Ar(L)exp(-A25/L)} . (28) 



Now 



1 A;2/23v^ 



4 V5 25 



and 



(29) 



C-^ {C [n{L) exp{-A^By/L)] } = N{L) exp{-A^By/L ] 



Using the convolution property, Eq. (|T3D, we arrive at 



{L-u)-'^(^U\L-ur/^ + ^{L-uy 



(30) 



N{u) eyiY>{-A^B^/u)du . 

(31) 



By examining Eq. (^Tj) we see that there are 3 integrals of the form 

I^ = c (L - uyN{u) exp{-A'^B^/u) du , 

JO 

with s = 1, s = 5/2 and s = 3. Writing N{u) from Eq. (^) as follows 



(32) 



N{u) = 1 /(Mu'/^ + Tu^/^exp{-A^B^) 



(33) 
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where M = 2A'k'^/'^ / \fB and T = tx/2A?B'^. Simple manipulations ( see Appendix C ) lead 
to the result 



(34) 



where the integral $„, which has some finite value in the large L limit, is given by 



lo T + Mx^/^exp{A'^Bx] 
Using these formulas we derive 



dx 



(35) 



Ji = 2L$o - 2$2 



'5/2 



-k 



2 / 2*^0^2 ^^2 



2$. 



L + ... 



(36) 



AB 2B 

Using Eqs. (|36| ) in Eq. (pT]), and keeping only the leading order terms in L, leads to the 
result 



G2(k,L) = 2L% - 2<!)2 - k^ 



2A^^5/2 ^ 3£o^2 _ ^2^^3/2 _ 3£2^ 

5 45 ^ 25 



(37) 



By substituting (H), (1^) and (0) in Eq. (0) we find 



{R') = -4- 



.1 (^^2^3/2 ^ ^ A^2^3/2 _ ^ (_M!i0^5/2 _ 3^^2 ^ ^2^^^3/2 ^ 



2B 



1 - cj(2L$o - 2<l>2) 



(38) 



-2a;$o^ + 2u;$2 + 1 

Here the integrals are all bounded and hence we can take the large L limit of the entire 
expression freely. Cancelling out the corresponding —2uj^L term we get 



^ ^ 5 2B 



3 , (2A^^2 



+ 



\ $0 8a;<l>o^ 
15 _ 20a^$2+|^_i 



(39) 



which is our final result. We see that in deriving this result, as no divergences are encoun- 
tered, there is no restriction on the value of w for the convergence of the series in the large 
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L limit. Comparing Eq. ( pQ]) with Eq. (|T]), we see that A = 1/2 and the coefficient C is 
negative. However, unless the calculation is taken to the next order, we cannot be certain 
about the sign of C as there are additional contributions for this term to 0(y{sY). We note, 
however, see section IV, that a negative C is consistent with our Monte Carlo simulations. 



IV. DISCUSSION 



The prediction A = 1/2, the main result of this paper, differs from all previous pre- 



dictions with the exception of that of Saleur [11| who rejected it for technical reasons. We 



summarise the various predictions for A in 2D in Table (|) . The numerical simulations are of 



particular interest. With the exception of a very few authors [0,^] (who however were not 
concerned with the correction-to-scaling terms), it appears few simulation studies have been 
done on the continuum, most previous work concentrating on 2D and 3D lattices. Simula- 
tions of chains in the continuum is more demanding computationally than for a lattice. We 



have addressed this problem with the use of a biased sampling Monte Carlo method p6|j27 
to create 2D self-avoiding walks on the continuum, and have applied a finite size scaling 
analysis to the resulting data. We have found A = 1/2 and a negative C coefficient. 



consistent with Eqn. (^). We note some inconsistencies with lattice simulations in view of 
the universality of A ||2^, indicating that the analysis of this data is strongly affected by 
assumptions regarding C in Eq. (|l|). These details will be published elsewhere |^ . 



Although the above calculation can also be performed starting from the solution to the 3D 



Edwards model |]T4[, resulting in a convergent series in A^, it is unlikely the correction term 
thus calculated would be reliable. We provide the following heuristic argument to support 
this. Since the 3D Edwards model prediction z/ = 3/5 is now known to be inexact |31[|, one 
can easily show from a simple dimensional analysis that for [r] ~ L^'^^ , where z/ = 3/((i + 2), 
our perturbation potential /g^ V{s) ds in Eq. (||) varies as L(^~'^)/'^'^+^)(aiL~^(^'^~^)/^— a;2-^^'^'^)- 
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Since v = 3/4 is exact in 2D, i.e., 5 = 0, it is possible, provided the amplitudes ai = a2 Q, 
for there to be a cancellation of divergences to all orders, resulting in controlled errors in the 
perturbation series. This, however, cannot be true in 3D, since 6 < 0, even if the amplitudes 
are equal. 

In conclusion, we have presented the results of a new perturbation method which starts 
from the Edwards 2D self-consistent solution. As the later already correctly predicts the 
exact exponent z/ = 3/4 for {R%), the resulting perturbation series is shown to be convergent 
in and free from uncontrolled terms. We find the leading correction-to-scaling exponent 
A = 1/2, and the next order coefficient C to be negative. These values are supported by an 



analysis of self-avoiding walks on the 2D continuum |^. It appears that A = 1/2 is also 



possible on the 2D lattice after a reassessment [|30| of current lattice data analysis methods, 
and that there is new physics in the 2D Edwards solution. 

TCC would like to thank Prof. Sir S. F. Edwards, Prof. D. Sherrington, Prof. R. 
Stinchcombe and Prof. M. Barma for helpful discussions. SRS acknowledges the support of 
an Australian Postgraduate Research Award. 

APPENDIX A: CALCULATION OF N{L) 

From Eq. (120|) we have 

J ^^(R, L)(fR = j^"" N{L) exp (-^(^ - AL^'^^ R dR dO 

B 



= 2TtN{L) Rexp [-j{R - AL^/^)^^ dR . (Al) 
Let X = i? - 



^This statement requires the equality of the exact amphtudes, ai and 02, as clearly the 'pertur- 



bative' (i.e. approximate) amplitudes cancel by virtue of the convergence of Eq. (|39D and its 3D 
equivalent. We are unfortunately not able to prove this at present. 
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2'kN(L) I ix + AL^/^)e--'' dx 



= 2ttN{L) 
For large x, the Erf(x) is given by 

Erf(x) = 1 



— e-^'^-^ + £lV^^5/4 Erf (Av^L^/^^ 
2B 2/B ^ ^ 



+ 



^/7^X ^^X^ 

Keeping the leading order terms, (|A1|) becomes 



J^^-A'^BVL , ^V^ r5/4 _ -A^Bv^ , -^^^^ 

2B Jb 2B ^ 2A^B^' 



2ttN{L) 

As noted in Eq. (]22|), this simply equals 1. Thus 

1 /("^l^^ + ^lVV^^^v/^' 



iV(L) 



(A2) 



(A3) 



(A4) 



(A5) 



APPENDIX B: CALCULATION OF G£;(K, L) 



From Eq. (|20|) we have 



2n roc 



JO 



ikR cos 9 ; 



N{L) exp --(i? - AL3/4)2 



27rA^(L) RJoikR) exp (-"^(^ - AL^/^^ dR 



To order fc^, Jo{kR) = 1 - 



4 ' 



thus 



Gi^fk,/.) =27rAr(L) / R(l — )exp( {R-AL^/y]dR 

Jo 4: \ L 

= 2itN{L) Rexp (^~{R - AL^/^Y^ dR 

-p|iv(L) R' exp [-jiR - dR . 



(Bl) 



(B2) 



The first part of integral simply equals 1, see Eq. p^), while the second k"^ dependent part 
equals 
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Let X = R — AL^^^ and the above integral becomes 

nm^A^B-/! ZLR3/2 V ^ 



-e-N{L) 



On expanding the error function, several leading terms cancel leaving 



,7r 



k'-N{L) 



B 253/2 3^454 



Thus 



TT 



Ge{KL) = i-r-iv(L) 



(B3) 



(B4) 



(B5) 



(B6) 



253/2 ■ 8^454' 

Substituting N{L) into the above and dropping exponentially small terms in the large L 
limit, we get 



2 2A7r3/2^ 



253/2 



5 



^9/4 



^ 4 



^2^3/2 + 

25 . 



(B7) 



APPENDIX C: CALCULATION OF THE INTEGRAL Is 



The general form for the required integral is given by 



L {L-uYexp{-A'^By^) 
^ Mm5/4 + Tu^/^expi-A^B^) ' 



By a series of simple substitutions we can simplify ( |Cll) to 



2cL' 



2cU 



T + Mx3/2exp(A25x) 
1 



dx 



dx — 



X 



(CI) 



/o T + Mx3/2exp(A25x) L 7o T + Mx3/2 exp(A25a;) 
2cL'% - 2csL'-^^2 + C(L""^$4) + ••• , 



(ia; + 



(C2) 
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where the integral $p is given by 

T + Mx3/2exp(yl25x) 

and is bounded for all L. 



" Jo T + Mx^/^e-KviA^Bx) ' ^^^^ 
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TABLES 



TABLE I. Predicted values for A in 2D 



Reference 



Method^ 



Predicted A 



Nienhuis @ 



Chaves and Riera [21 



Saleur [|ll[ 
Havlin and Ben-Avraham 
Rapaport Q 
Wang m 
Lyklema and Kremer ^ 
Djordjevic ei a/ ||5| 
Privman ^ 
Ishinabe @ 
Lam [23 
This work 



Shannon et al [30 



Hexagonal lattice mapping 
Cell renormalization group 
Conformal invariance 
MC= 

Exact triangular*' 

MC^ 
Exact triangular^ 
Exact triangular^ 
Exact square' 

MCJ 

Improved perturbation 
MC continuum'^ 



^MC = Monte Carlo, ESS = Finite size scaling 
''This result was rejected. 

*^Square lattice data using extrapolation of series analysis 
"^Square and triangular lattice data using linear regression 
^'Fitting to a series analysis 

^Square lattice data using extrapolation and fitting to a series analysis 
^Extrapolation and fitting to a series analysis 
'^Intersection of a ESS series analysis 
'Intersection of a ESS series analysis 

j Square and triangular lattice using integration of a ESS series 
'^Fitting of a ESS series analysis 



3/2 
1 

11/16, 1/2^^ 
1.2 ± 0.1 
1 

0.85 ± 0.05 
0.84 

0.66 ± 0.07 
0.65 ± 0.08 
0.65 ± 0.05 
0.6 
1/2 

0.5 ± 0.05 
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